The Crab Nebula's Moving Wisps in Radio 
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ABSTRACT 

We present three high resolution radio images of the Crab nebula, taken in 1998.6, 1998.8 
and 2000.1 with the VLA. These are the best radio images of the Crab to date. We show that, 
near the pulsar, there are significant changes between our three observing epochs. These changes 
have an elliptical geometry very similar to that of the optical wisps. One radio wisp in particular 
can be unambiguously identified between two of our observing epochs, and moves outward with 
an apparent velocity of 0.24c. The similarity in both morphology and behavior of the present 
radio wisps to the optical wisps suggests that they are associated. This implies that the radio 
wisps, like the optical ones, are likely manifestations of the shock in the Crab pulsar's wind. 
This suggests that the radio emitting electrons are accelerated in the same region as the ones 
responsible for the optical to X-ray emission, contrary to most current models. 

Subject headings: ISM: Individual (Crab Nebula) — Radio Continuum: ISM — supernova remnants 



1. Introduction 

The Crab Nebula is one of the most important 
objects in astrophysics. It is a beautiful example 
of a pulsar powered nebula, and is also the most 
accessible pulsar nebula. Each generation of ob- 
servations has revealed new intricacies to its struc- 
ture. Recent results have highlighted the remark- 
able nature of the region surrounding the Crab 
Nebula's pulsar, where the pulsar wind is injected 
into the nebula. A remarkable image of this region 
from the Chandra X-ray observatory (Weisskopf 
et al. 2000) shows a complex geometry consisting 
of a tilted toroidal structure with a jet along the 
symmetry axis of the torus. A recent sequence of 
images from the Hubble Space Telescope (HST) by 
Hester et al. (1996) revealed astonishing details in 
this region. In particular, there is a series of ellip- 
tical ripples in this region, which are usually called 



"wisps" . 

The wisp region has long been known to be 
rapidly variable (Lampland 1921; Scargle 1969), 
but only recently has the nature of the variability 
been conclusively established: Hester et al. (1996; 
see also Hester 1998) and Tanvir, Thomson, & 
Tsikarishvili (1997) have demonstrated that the 
wisps are variable on time-scales as short as a few 
days and that they are moving outwards with ap- 
parent speeds of up to 0.7c. Early Chandra results 
indicate that the X-ray structure is also variable, 
on timescales of ^ 10 days. 

At the center of the Crab Nebula, and of the 
wisp region, the rapidly spinning Crab pulsar is 
slowing down, and thus losing energy at the rate of 
5 X 10"^^ erg s~^. This bulk of this energy emerges 
from the pulsar as a highly coUimated wind of rel- 
ativistic particles and magnetic field. This outflow 
is randomized at a shock, where it comes into equi- 



librium with its surroundings (Rccs & Gunn 1974) 
and then flows into the body of the nebula, where 
it gives rise to the synchrotron emission we ob- 
serve from the body of the nebula. However, the 
process by which this transfer of energy happens is 
not yet well understood (see e.g. Begelman 1999; 
Arons 1998; Hester 1998; Chcdia et al. 1997). 

Since almost all the pulsar's spindown energy is 
accounted for in the body of the nebula, we know 
that the transfer of this energy is quite efficient. 
Only a very small fraction of the energy seems 
to be radiated directly from the region where the 
wind is randomized: the wisps are now thought to 
be associated with the shock which terminates the 
highly collimated outflow from the pulsar and ran- 
domizes it before it enters the body of the nebula. 
These features are our only observational window 
into the process by which the pulsar outflow is con- 
verted into the rclativistic gas and magnetic fleld 
which fill the body of the nebula. 

A feature with an arc-like geometry similar to 
that of the optical wisps was observed in the radio 
spectral index by Bietenholz & Kronberg (1992), 
who also observed changes in the radio emission 
from the center of the nebula over a period of five 
years. Our goal with these observations is to inves- 
tigate the nature of the radio synchrotron emission 
from the wisp region, and to compare it to the syn- 
chrotron emission from the optical wisps. In par- 
ticular, we will investigate whether the radio wisps 
are also rapidly variable, since the extant observa- 
tions of them show only that they are variable on 
a 5 year timescale. 

2. Observations and Data Reduction 

We observed the Crab Nebula using the NRAO 
VLA^ in the 5 GHz band. Table 1 shows details of 
our three observing runs, all of which were done 
using the VLA in the B configuration. Due to 
scheduling reasons, our observing run in February 
2000 was split into two segments, which we com- 
bined, and we refer to this combined data set as 
our 11 February 2000 epoch. At each of these ses- 
sions, we observed at a spaced pair of frequencies 
(within the 5 GHz band) in order to increase our 
u-v coverage, since good u-v coverage is critical 

^The NRAO Very Large Array is a facility of the National 
Science Foundation operated under cooperative agreement 
by Associated Universities, Inc. 



when imaging an object as extended as the Crab. 
In addition, we chose different pairs of frequencies 
for the 10 Feb. and the 12 Feb. 2000 observing 
runs, since these two runs had roughly the same 
range in hour angle. 

Our primary goal was to investigate any rapid 
variations in the nebular structure. Since a speed 
of c represents a proper motion of only ~ 3" per 
month, we do not expect the large scale struc- 
ture of the nebula to change very much on short 
timescales. We do know that the radio syn- 
chrotron nebula is expanding at 0.13% year~^ (Bi- 
etenholz et al. 1991). The effect of this expansion 
is both predictable, since the rate is known, and 
is also quite small over the timespan of our obser- 
vations: the most rapid motion due to the overall 
expansion rate is < 0'.'3 yr~^. 

The fact that the large scale structure of the 
nebula is not expected to change over our times- 
pan is crucial to our strategy: an interferometer 
like the VLA measures the spatial Fourier trans- 
form of the sky brightness, and is sensitive only to 
a range of spatial frequencies determined by the 
length of the interferometer baselines and the ob- 
serving frequency. For example, in the B configu- 
ration at 5 GHz, the VLA is sensitive only to struc- 
ture on spatial scales between ~ 1'.'3 and ~ 50". 
The Crab Nebula, however, exhibits structure at 
spatial scales from ^ 6' down to < 1". Typically 
in such cases, one would obtain the larger scale 
structure from VLA observations in more compact 
array configurations. However, since VLA configu- 
ration changes occur only once every four months, 
that approach is not compatible with our goal of 
obtaining time-resolved observations. 

Our solution to this problem is to use maxi- 
mum entropy deconvolution (AIPS task VTESS; 
see Cornwcll & Evans 1985, Cornwell 1988), which 
allows us to recover the large scale structure by 
supplying a low resolution support (default) in the 
deconvolution process. The deconvolved image is 
then biased to be as close to the support as is al- 
lowed by the data. We chose this support to be the 
same for all three of our epochs^ which serves to 
make any differences between the epochs be only 
such as are demanded by the data. 

As our support, we used an image of the Crab 

^Wc did scale the support according to the known expansion 
of 0.13% year~i but the effect of this is minimal. 
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made using data taken in 1987, at 5 GHz, using 
the B, C and D configurations of the VLA (Bieten- 
holz & Kronberg 1990; 1991). We convolved this 
image with a round Gaussian of full- width at half- 
maximum (FWHM) 20", and scaled it to account 
for the general expansion of the radio nebula since 
1987. Because wc have only B array data at the 
current epochs of interest, the structure in our im- 
ages on the largest spatial scales, i.e. those > 1', 
is derived predominately from the support image, 
and we thus have very little information on any 
changes that might occur at these spatial scales. 
As we have argued above, however, there is little 
reason to expect rapid variations on these scales. 
On spatial scales smaller than ~ 30" on the other 
hand, our 21-v coverage is excellent, and wc believe 
that our images reliably indicate the structure of 
the nebula, and in particular, changes therein from 
one epoch to another. 

3. Results 

An image of the Crab Nebula on 11 Febru- 
ary 2000 is shown in Figure 1. The restor- 
ing beam size was 1'.'4 (this restoring beam size 
was conservatively chosen as a common size for 
all three epochs). The background rms level 
was 76 /LtJy beam~^ before applying the pri- 
mary beam correction, and the peak flux was 
46.6 mJy beam~^. This is the highest dynamic 
range and resolution radio image yet produced of 
the Crab Nebula. We note that it is superior to 
the three and four-configuration VLA images of 
Bietenholz & Kronberg (1990; 1991) — on which 
it is partly based through the use of those images 
as a support — because of our presently increased 
u-v coverage at the higher spatial frequencies. 

The images taken 1998 August and October 
were very similar, and so we do not reproduce 
them in their entirety. They had background rms 
levels of 100 and 82 /iJy bcam^^ respectively. In- 
stead, we present a difference image: Figure 2 
shows the difference between the images of Oc- 
tober and August 1998. 

On this difference image a prominent series of 
elliptical ripples are visible near the center of the 
nebula. The brightness of these ripples is several 
times that of the background noise level. Such 
ripples arc only visible near the pulsar position: 
Figure 2 shows that, in the rest of the nebula. 



any difference features are of lower amplitude and 
larger spatial scale, and likely due to noise, missing 
short-spacings, and errors in the deconvolution. 
The elliptical ripples show a remarkable similar- 
ity in geometry with the elliptical wisps seen in 
the optical (Hester et al. 1995; Hester et al. 1996) 
and the torus visible in the X-ray (Aschcnbach & 
Brinkman 1975; Weisskopf et al. 2000). Figure 3 
shows a region near the pulsar of the radio differ- 
ence image between 1998 October and 2000 Febru- 
ary, and a corresponding difference image between 
two HST images taken 1.5 months apart and con- 
volved to the same resolution (note that the im- 
ages in this figure have been rotated by —47°). In 
particular, the implied inclination angle (assuming 
that the features are intrinsically circular) is sim- 
ilar to that in both the optical and the X-ray. If 
the ripples were an artefact, such a coincidence in 
geometry with the real features observed at other 
wavelengths would be quite unlikely. While mov- 
ing features seem to exist over roughly the same 
region in both the optical and radio images, the 
largest features in the radio are further from the 
pulsar than the largest optical features are. 

In order to make the motions from one epoch to 
the next better visible in reproductions, we show 
in Figure 4 high-pass filtered views of the center 
of each of our three images. We used a Gaus- 
sian high-pass filter with a FWHM of 25", which 
serves to accentuate the relatively faint wisp-like 
mobile features. We will call the most mobile fea- 
ture wisp a, and we indicate its position in 1998 
August by a black line on all three images in or- 
der to show the displacement in the subsequent 
epochs. 

These motions arc real: they occur on spatial 
scales of a few arc-seconds, exactly those to which 
we are sensitive and for which we have excellent 
u-v coverage. In order to further ascertain that 
these motions were not some form of artefact due 
to the deconvolution process, we performed several 
tests. Firstly, we deconvolved the 1998 August im- 
age using the October image as a support, instead 
of our usual low-resolution image. This further' bi- 
ases the images to be as similar as is allowed by the 
data. The ripples remained essentially unchanged. 
Secondly, expanding or contracting the default im- 
age by 2% also left the ripples unchanged, but re- 
sulted in a poorer image. Lastly, we self-calibrated 
the 1998 August u-v data both in amplitude and 
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Table 1 
Observing Run Details 



Date 


Length 


Frequencies 


Midpoint 




(hrs) 


(MHz) 


(days since 1998.0) 


1998 Aug 9 


9.6 


4615, 4885 


221.7 


1998 Oct 13 


10.1 


4615, 4885 


286.5 


2000 Feb 10 


3.7 


4749, 4996 


772.1 


2000 Feb 12 


5.7 


4615, 4885 


774.0 


(combined Feb 2000) 


9.4 


4615, 4749, 4885, 4996 


773. P 



''The midpoint for the combined February 2000 runs was taken to be the weighted 
mean of the midpoints of the two individual days. 



phase, using the October image as a model, and 
vice versa, and then re-imaged. Again, the rip- 
ples are essentially unchanged, implying that they 

cannot be ascribed to any calibration anomalies. 

Finally, in order measure the speeds of fea- 
tures, we plotted profiles through the high-pass 
filtered images. Figure 5 shows profiles for our 
three epochs, drawn through the pulsar at posi- 
tion angles, p. a., of —11°, chosen to pass through 
the most mobile wisp, and at —90°, chosen so as 
to pass through a mobile feature on the opposite 
side of the pulsar. There are numerous changes in 
addition to the motion of wisp a indicated in Fig. 4 
above. However, with our very coarse sampling in 
time, the identification of features from one epoch 
to the next is difficult. We note that the pulsar, 
which is known to be quite variable but has a mean 
flux density of 1.4 mJy at 4.8 GHz (Lorimer et al. 
1995), is visible as a variable peak at the origin in 
the profiles. 

The identification is wisp a between 1998 Au- 
gust and October is unambiguous. It moves with 
an apparent velocity of 0.24 ± 0.06c. We indicate 
this motion with a solid line on Fig. 5a. By 2000 
February 11, the identification is no longer com- 
pletely unambiguous — we mark the motion de- 
rived from the likeliest identification in Fig. 5a 
with a dashed line. This identification implies 
an average apparent speed between 1998 August 
and 2000 February of 0.28 ± 0.01c, which is con- 
sistent with the speed between 1998 August and 
October. On the opposite side of the pulsar, the 



most prominent moving feature is seen on the 
profile at p. a. = —90°. We call this feature the 
counter-wisp, and mark its motion in Fig. 56. The 
speed between 2000 August and October is not 
reliably determinable, but the apparent speed be- 
tween 1998 (taking the average of August and Oc- 
tober) and 2000 February is 0.08 ± 0.02c, and this 
motion marked in Fig. bh. We note that because 
of the large time interval and the complexity of 
the changes between our 1998 and 2000 epochs, 
the identification of features over this time-span is 
necessarily somewhat more tentative. 

Assuming an inclination angle of 30°, we can 
calculate the true speed of the features above: 
wisp a has a true speed of 0.38 ± 0.01c, (assum- 
ing that the above identification in Feb. 2000 is 
correct; see Fig. 5a). The counter- wisp has a true 
speed of 0.19 ± 0.05c. 

4. Discussion 

We have shown that there are rapidly moving 
features, with velocities < c/3, visible in the radio 
emission from the Crab Nebula. These features 
have an arc-like geometry similar to that of the 
optical wisps. Bietenholz & Kronberg (1992) had 
already reported the existence of similarly shaped 
features near the center of the nebula, apparent in 
the radio spectral index. They noted the possibil- 
ity that these features might in fact be temporal 
changes rather than true spectral index features, 
because the two images from which they derived 
the spectral index were effectively taken several 
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months apart"^. In light of the present results it 
seems likely that the features reported in Bieten- 
holz & Kronberg (1992) were in fact the result of 
temporal changes, rather than true spectral index 
anomalies. We cannot of course exclude the pres- 
ence of true spectral index variations, but the data 
so far at hand do not demand them. 

We propose that these rapidly variable radio 
wisps are associated with the similar features seen 
in the optical and the torus visible in the X- 
ray. The precise geometry of the radio features 
is harder to discern than that of the optical wisps 
due to the presence of a relatively much brighter 
complex background provided by the rest of the 
nebula. Nonetheless our difference images give a 
distinct picture of the temporal changes in these 
features, and clearly show an elliptical geometry. 
In fact, the ellipticity is the same as that of both 
the optical wisps and the X-ray torus, suggesting 
that, like the optical wisps and the X-ray torus, 
the radio wisps are circular features seen at the 
same inclination angle of ~ 30° . 

Like the optical wisps, the radio wisps seem to 
be moving outward (although this is not clearly es- 
tablished by our observations, since we have only 
three epochs and sub-optimal time sampling, the 
unambiguously identified wisp a is moving out- 
wards between August and October 1998). Fur- 
thermore, we measure speeds of < c/3 for the 
features which wc can identify from one epoch to 
the next. This is in fact the expected speed for 
the post-shock flow (Blandford & Rees 1974) and 
consistent with the speeds measured in the opti- 
cal (Hester et al. 1996; Hester 1998; Tanvir et al. 
1997) and inferred from the X-ray (Felling et al. 
1987). These arguments strongly suggest that the 
presently observed radio wisps are associated with 
the optical wisps and the X-ray torus. 

In order to discuss the implications of these ob- 
servations, we will briefly describe our theoreti- 
cal understanding of the Crab Nebula. As men- 
tioned in § 1, the visible synchrotron nebula is 
powered by the spindown energy of pulsar. This 
energy emerges from the pulsar in the form of 
a collimated, highly relativistic wind, which con- 

^The spectral index was between 5 GHz and 1.4 GHz. The 
5 GHz B-ajrray data were taJsen on 1987 November 25, while 
the resolution-matched 1.4 GHz A-array data were taken 
on 1987 July 7. 



sists of some combination of particles and mag- 
netic field. The coupling to the outer boundary 
condition, namely the outside of the observable 
synchrotron nebula which moves with f <C c, re- 
quires a shock in the wind. It is at this shock that 
the particle velocities are randomized, and it is the 
post-shock wind which then supplies the relativis- 
tic particles and the magnetic field responsible for 
the nebular synchrotron emission. The wisps are 
thought to be manifestations of this termination 
shock, since pressure balance arguments place the 
location of the shock roughly at the wisp location 
(Rees & Gunn 1974). Similar features have been 
seen near other pulsars (Bietenholz, Frail, & Han- 
kins 1991; Frail & Mofi'ett 1993), and so they may 
be a common manifestation of the interaction of a 
pulsar wind with its surroundings. In particular, 
the wisps are thought to be associated with the 
shock in the equatorial sector of the pulsar wind, 
since both the optical and especially the X-ray ob- 
servations (Weisskopf et al. 2000) clearly show a 
toroidal geometry. (We note that polar features 
are also seen in both the optical and the X-ray, 
but as yet we cannot conclusively identify any po- 
lar features in the radio). 

The precise mechanism of wisp formation, how- 
ever, is not yet well understood (Bietenholz & 
Kronberg 1992; Chedia et al. 1997; Hester 1998; 
Arons 1998; Begelman 1999). Also, while the ra- 
dio spectral index is very uniform over the neb- 
ula (Bietenholz et al. 1997), suggesting a single 
and presumably central source for all the radio 
emitting electrons, their origin in the shock near 
the pulsar is not yet understood (see e.g., Atoyan 
1999; Arons 1998). At the termination shock of 
the wind, one expects that the mean Lorentz fac- 
tor, 7, of the post-shock particles to be roughly 
equal to that in the wind, which is widely thought 
to be ~ 10*^ (Rees & Gunn 1974; Kennel fc Coro- 
niti 1984; Melatos & Melrose 1996; Arons, 1998). 
However, this would produce essentially no radio- 
emitting, i.e. 7 ~ 10^ electrons. 

For example, the model of Arons (1998; see also 
Gallant & Arons 1994) can account for wisps and 
also the observed nebular spectrum from the op- 
tical to the X-ray, but it cannot account for the 
radio-emitting electrons. Arons suggests that a 
slower, high-latitude wind might produce the ra- 
dio emitting electrons. In the model of Atoyan 
(1999), which can account for the entire nebular 
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spectrum, the radio emitting electrons are of his- 
torical origin, having been emitted when the pul- 
sar was younger and spinning much more rapidly. 
However, in neither of these cases would radio- 
emitting electrons be accelerated in the current 
wisp region. Thus the detection of wisps in the ra- 
dio implies that the acceleration of particles in the 
equatorial shock termination region produces low 
energy radio-emitting electrons as well as those 
at higher energies responsible for the optical and 
higher synchrotron emission. 

In summary, we present new radio images of the 
Crab Nebula, with higher resolution and dynamic 
range than those previously available. We show 
that there are significant changes occurring over 
periods of two months and one year. An arc-like 
feature, or wisp, is seen to be moving outwards, 
and images of the differences between successive 
pairs of our three observing epochs clearly show 
numerous features which have an elliptical geom- 
etry. We propose that these radio wisps are asso- 
ciated with the optical wisps, which suggests that 
acceleration mechanism producing radio-emitting 
electrons is the same as, or at least co-spatial to, 
that producing the more energetic electrons re- 
sponsible for the nebular optical and X-ray syn- 
chrotron emission. Future, more highly time- 
resolved observations will hopefully allow us to es- 
tablish the exact nature and speed of the motions. 

Research at York University was partly sup- 
ported by NSERC. We thank Barry Clark for his 
patience with the scheduling difficulties involved 
in these observations. 
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Fig. 1. — The image of the Crab at 5 GHz on 2000 February 11 (from data taken on Feb. 10 and 12), 
after primary beam correction. The FWHM size of the restoring beam was 1'.'4. The peak flux density is 
46.6 mJy bearn"^, and the exterior rms was 76 /zJy beam~^ before primary beam correction. The contours 
are drawn at 0.75, 2, 4 and 8% of the peak, and the greyscale is in mJy beam"""^. Maximum entropy 
deconvolution was used with a support to recover the low spatial frequency structure (see text, § 2) for 
details). 
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BIGHT ASCENSION |[B1 950) 

Fig. 2. - The differences between our 1998 October 13 and 1998 August 9 images in greyscale, with the 
grayscale ranging from +5 to —5 mJy beam~^. The FWHM size of the restoring beam was l'/4. For 
reference, we also show the the 2, 5, and 10% contours from the 2000 February 11 image, and indicate the 
pulsar position by a cross. The box indicates the region of the high-pass filtered images shown below in 
Fig. 4 
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Fig. 3. — The difference features in the optical and radio. The images have been rotated by —47°. The cross 
is drawn in the N-S E-W direction, and shows location of the pulsar. In a) wc show the radio difference 
image between 2000 February 11 and 1998 October 13 (the greyscale ranges from +5 to —5 mJy beam~^, 
and the convolving beam size was 1'.'4 FWHM). In b) we show the difference between the HST images from 
2000 October 23 and September 9, convolved also to 1'.'4 FWHM resolution. 
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Fig. 4. — High-pass filtered images of the region near the pulsar, showing the mobile wisps. The images 
have been high-pass filtered with a Gaussian of FWHM 25". The restoring beam again has a FWHM size of 
1'.'4, shown at lower left. The position of wisp a on 1998 August 9 is indicated on all three images to show 
its displacement. In a) we show the image from 1998 August 9, in b) that from 1998 October 13, and in c), 
that from 2000 February 11. 
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Fig. 5. — Profiles through the center of the Crab at three epochs. The profiles were made from the high- 
pass filtered images shown in Fig. 4. In each case the bottom profile is that from 2000 February 11. The 
middle and the upper profiles arc those from 1998 October 13 and August 9, respectively, and have been 
artificially shifted upwards by an amount proportional to the the time difference to the 2000 February epoch. 
The profiles are drawn through the pulsar position (5^ 31" 31!423, 21° 58' 54'.'26), which is indicated by the 
dotted vertical line. In a) wc show profiles drawn at a p. a. of —11°, with a +ivc displacement being to the 
NW. The motion of wisp a is indicated, with the solid line indicating the secure identification of this feature 
between 1998 August and October, and the dashed line showing the likeliest identification in 2000 February. 
In b) we show profiles in R.A. or at a p. a. of 90°. The motion of the counter- wisp is indicated. 
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